mdot=8.5:0.01:100;
b1=73./mdot.^2;
b2=(1-1./(sqrt(1+2.*73./mdot.^2)));
hold off
plot(mdot,b1)
hold on
plot(mdot,b2)
hold off
% 17 beaming factor 16 m_tr/m_edd
% 3 Lx
minLx=37;
[all_rb_1,all_wind_1]=DELallrbwind(all_rb,all_wind);

beaming_ori=all_rb_1(:,17);
beamingmodel1=getBeaming( 1,all_rb_1(:,16));
beamingmodel2=getBeaming( 2,all_rb_1(:,16));

subplot 312
plot_beaming_gedian(all_rb_1,0);
title('origin')
all_rb_1(:,17)=beamingmodel1;


subplot 313
plot_beaming_gedian(all_rb_1,0);
title('model 1 ')

all_rb_1(:,17)=beamingmodel2;
subplot 311
plot_beaming_gedian(all_rb_1,0);
title('model 2 ')